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ABSTRACT 

The spectrum of the Io plasma torus in the range 905-1187 A was recorded at 0.26 A 
resolution by the Far Ultraviolet Spectroscopic Explorer (FUSE) on 2001 January 14. 
Five orbits of data were obtained with the west ansa of the torus centered and tracked in 
the 30" x 30" apertures of the FUSE spectrographs for a total observation time of 9740 
seconds. This region of the spectrum is dominated by transitions of S II, S III and S IV, 
whose multiplet structure is nearly completely resolved. We confirm our earlier detection 
of emission from resonance multiplets of CI III and CI II and derive an abundance of 
Cl +2 of 2.1% relative to S 2+ , leading to an overall chlorine ion abundance in the torus 
of slightly less than 1%. A number of features near 990 A remain unidentified, and C III 
A977 is detected in two independent channels at the 3-a level. The inferred relative ion 
abundance of C 2+ relative to S 2+ is 3.7 x 10 -4 . We also present spectra at 0.085 A 
resolution taken on 2001 October 19 and 21 with the 4" x 20" aperture. In these 
spectra the observed lines are resolved and their widths correspond to ion temperatures 
of 60-70 eV for all three sulfur ions. 

Subject headings: planets and satellites: individual (Io, Jupiter) — ultraviolet: solar 
system 



1. INTRODUCTION 

In a previous Letter (Feldman et al. 2001a, hereinafter Paper I) we reported on early spectro- 
scopic observations of the Io plasma torus with the Far Ultraviolet Spectroscopic Explorer (FUSE). 
Those observations were made in January 2000 before FUSE acquired the capability to track mov- 
ing solar system targets yet yielded new information about the content of the torus and in particular 
allowed the determination of the relative abundance of chlorine ions from the first detection of the 



x and Department of Earth and Planetary Sciences, The Johns Hopkins University 

2 present address: Space Department, Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins 
Road, Laurel, MD 20723-6099 



- 2 - 



ultraviolet resonance transitions of both CI II and CI III. The source of chlorine on Io is currently 
of interest following the detection of NaCl in Io's atmosphere (Lellouch et al. 2003) and recent de- 
tailed modelling of volcanic chemistry (Moses et al. 2002). As part of the FUSE team solar system 
project, the torus was observed again, this time with full tracking capability, on 2001 January 14, 
during a Jovian campaign timed to coincide with the near closest approach of Cassini to Jupiter. 
The spectra obtained were a significant improvement over the earlier data. A clear enhancement of 
signal-to-noise ratio by a factor of four enabled confirmation of the identification of the weak CI II 
emission reported earlier, the detection of many weak S II and S III emissions including several from 
very highly excited states, and the spectrum between 905 and 990 A that was not obtained from the 
earlier data due to thermal misalignment of the SiC channels. Additional, higher resolution spectra 
of the strongest emissions were obtained during two series of limb scans of Jupiter in October 2001. 

In this paper we present a more detailed analysis of the spectrum and in particular examine 
the torus abundances that can be inferred from the presence or absence of emissions of ions of C, 
N, Si and P, all of which have strong transitions in the FUSE spectral range. A number of weak 
spectral features remain unidentified and there are a few cases where the spectral assignment of the 
observed lines are uncertain. The higher resolution spectra show line widths that are greater than 
the instrumental widths which allows us to directly determine the ion temperature in the warm 
region of the torus where the ultraviolet emissions are produced. A comparison with the earlier data 
shows a significant change in relative ionization stages of sulfur in the year between observations. 
Temporal variations on scales of days to months have been reported from extended Cassini UVIS 
observations made around the time of Jupiter fly-by (Steffi et al. 2003). The abundance of chlorine 
ions relative to those of sulfur also appears to have changed. 



2. OBSERVATIONS 

Observations of the west ansa of the torus were made on 2001 January 14, beginning at 05:56:40 
UT. A total exposure of 9740 seconds was obtained by accumulating spectra over 5 contiguous 
orbits. The center of the 30" x 30" aperture (LWRS) was placed at the position of the ansa at 
the beginning of each exposure and both the motion of Jupiter on the sky and the wobble of the 
torus were tracked during each individual exposure. Tracking was verified from the constancy of 
the count rates in the principal ion emissions over each ~2000 s exposure derived from the time- 
tagged data files. There is no spatial information about the distribution of the emission within the 
aperture. 

At the time of the observation, Jupiter's heliocentric distance was 5.05 AU and its geocentric 
distance was 4.40 AU. Details of the FUSE instrumentation have been given by Moos et al. (2000) 
and Sahnow et al. (2000). Basically, FUSE consists of four separate telescope/spectrograph as- 
semblies, two employing lithium fluoride coatings (LiF) and two using silicon carbide (SiC). Each 
of four separate detectors (denoted la, lb, 2a, and 2b) records two spectra, one each from a LiF 
and a SiC channel spanning a given wavelength interval, giving a total of eight separate spectra. 
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Due to the extended but finite size of the torus emitting region within the aperture, the effective 
spectral resolution is 0.26 A. The data presented here were reprocessed with the FUSE pipeline, 
calfuse version 2.2.1, all five orbits of data were co-added together, and the extracted fluxes were 
converted to average brightness (in rayleighs) in the 30" x 30" aperture. The wavelengths were 
reduced to laboratory values by first taking out a heliocentric velocity correction applied by the 
pipeline and then removing the doppler shifts due to the Earth's velocity relative to Jupiter and 
the plasma torus co-rotation of 75 km s _1 . The agreement with laboratory wavelengths in all cases 
is better than 0.05 A, and testifies to the accuracy of the wavelength calibration achieved in this 
version of the pipeline. 

Drift scans of the eastern limb of Jupiter were made on 2001 October 19 and 21 using the 
FUSE MDRS (4" x 20") aperture. A single scan was made over each of five consecutive ~3000 s 
orbits on both days. Each scan began with the slit, with the long dimension aligned nearly north- 
south, placed on a chord centered on the Jovian equator near the eastern limb. The spacecraft 
then executed a slew so that the apparent motion of the slit on the sky was towards the east, 
perpendicular to the long dimension of the slit, at a rate of ~0.23 arc-seconds per minute. At the 
end of each scan the slit was completely off the disk of the planet. For these observations, Jupiter's 
heliocentric distance was 5.14 AU and its geocentric distance was 4.86 and 4.83 AU, respectively, 
on the two dates. 

The goal of these observations was to record the spectrum of the disk, primarily resonance 
scattering by atomic hydrogen and electron impact of H2. This approach was chosen to minimize 
the possibility of contamination by the much stronger polar aurorae and to take advantage of 
possible limb brightening of the H I emissions. In this respect, the foreground torus emissions were 
contaminants as they had been in the lower resolution (~3 A) observations of Jovian airglow by 
the Hopkins Ultraviolet Telescope (Feldman et al. 1993). At the higher resolution of FUSE they 
are easy to separate from the disk emissions. 

Due to thermal misalignments between channels only the the LiFl channel appears to show 
the expected level of disk emission. This is the channel that is coupled to the Fine Error Sensor 
used for guiding and tracking. The SiCl channel also shows disk emission but at considerably 
lower signal/noise as the effective area of this channel is a factor of 4 lower than that for LiFl. The 
LiF2 and SiC2 channels both appear to have been off the disk for the entire observation, showing 
only torus emissions. For the LiFl channel the time tagged data clearly show the disk emissions 
decreasing with time while the torus emissions increase as more and more of the background 
emission becomes visible in the aperture. It is also possible to follow the brightness variation of 
the torus emissions over successive orbits as the centrifugal equator of the torus wobbles through 
Jupiter's equatorial plane at the limb. Because the emissions are about a factor of two weaker 
than at the ansae of the torus and because of the lower throughput of the MDRS aperture the 
torus spectra presented below are derived by co-adding all of the five orbits on each day for a total 
integration time of 31,333 seconds. 
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3. SPECTRA 

3.1. LWRS Spectra 

An overview spectrum from 995 to 1105 A obtained from the January 2000 observations was 
given in Paper I. The spectra from January 2001 are essentially similar except for the factor of 
four increase in signal-to-noise ratio. This is particularly important for confirming the CI II and 
CI III multiplets which are shown in Figure 1. The chlorine emissions, both of which are principal 
resonance transitions, occur close to the wavelengths of resonance transitions of S III and S IV. 
The components of the CI III triplet appear to be in the ratio 1:2:3 expected for the optically 
thin branching ratios even though the brightness of the 1015 A line is somewhat uncertain due to 
blending with the strong S III lines nearby. Similarly, five of the six lines of the CI II multiplet 
are clearly identified, again with the expected branching ratios. A list of all of the identified sulfur 
and chlorine ion emissions is given in Tables 1 and 2. Note that in two instances the atomic 
configurations are either not given or are uncertain. 

The tables also give the identifications and brightnesses of a number of lines not identified in 
the January 2000 spectra (Paper I), particularly of S II. This is not only due to the enhanced S/N 
but also to the absence of thermal misalignment between the LiF and SiC channels that occurred 
in the earlier observation. A segment of the SiC lb spectrum, which extends to the FUSE short 
wavelength limit at 905 A and includes the second resonance multiplet of S II, is shown in Figure 2. 
Figure 3 shows a region of the spectrum from the LiF2a channel containing highly excited multiplets 
of S III and S IV that are extremely useful as temperature diagnostics (see Section 4.2). 

3.1.1. Unidentified Features 

The high sensitivity and spectral resolution of this data set make possible the search for and 
identification of many other ionic species that have resonance transitions in the FUSE spectral 
range. A first examination of the data indicates that almost all of the observed emissions can be 
identified with known and published transitions of sulfur and chlorine. There are a few lines whose 
assignment is uncertain such as the S III line at 1016.85 A seen in the top panel of Figure 1 and 
listed in Table 2. Kelly (1987) gives a wavelength that matches the observed line very well but does 
not assign the transition. A weaker line seen in the same figure at ~1012.0 A is also listed by Kelly 
as a S III line at 1011.91 A, again without an assignment. The other ambiguity is the line at 1088.8 
A identified by Kelly as S IV. However, as we noted in our earlier paper, Kelly assigns this line to 
the same transition as the 1098.9 and 1099.5 A lines of S IV given in the NIST Atomic Spectra 
Database 3 . Of these latter two lines, Kelly lists only the 1099.5 A line, and does not identify the 
transition, so his assignment of the 1088.8 A line is suspect. The NIST Database does list a S IV 
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transition at 1088.214 A between high-lying 4 D° and 4 F states, the upper one requiring ~39 eV of 
excitation energy. However, the wavelength differs by 0.6 A from our measured wavelength and a 
companion line at 1091.18 A is absent, so this identification is unlikely. 

The only truly unidentified features in the spectrum occur near 990 A. Fi gure 4 shows the 
spectral region from 970 to 1000 A from three separate channels including the very short wavelength 
range of the LiFla channel. Five distinct lines are seen in the SiC2a spectrum just shortward of 
the S II line at 996.0 A, and the other spectra, where the wavelengths overlap, show the same 
features although at poorer S/N. These lines are at 991.5, 992.9, 993.5, 994.5, and 995.1 A and 
have brightnesses in the range of 0.3-0.6 rayleighs. It is tempting to assign the feature at 991.5 A 
to N III A991.56 but this is questionable due to the absence of the companion N III line at 989.79 
A which should be of comparable brightness. Kelly lists an unassigned S II transition at 991.73 A, 
but the wavelength agreement is not quite as good. No obvious single candidate multiplet for all 
of the lines has been found in the available spectroscopic databases. Of the other sulfur ions, there 
is a high-lying transition of S V at 994.46 A but this is considered unlikely since it requires 56 eV 
excitation energy and S 4+ is present only at about 2% relative to S 2+ (Steffi et al. 2003). There is 
a Si III triplet with lines at 993.5 and 994.8 A, close to the wavelengths of two observed lines, but 
the strongest member of the triplet should be at 997.4 A, which is missing. There should also be 
a stronger Si III multiplet at 1108-1113 A. 



3.1.2. Other Expected Features 

Both singly- and doubly- ionized carbon have strong transitions in the FUSE spectral range 
and these were searched for in the January 2001 data. The C III line at 977.02 A appears to be 
present in both SiC channel spectra (Figure 4) at the 3-cr level. Combining the two spectra gives 
a brightness of 0.19 ± 0.034 rayleighs. The C II doublet at 1036.34 and 1037.02 A is not detected 
with a 1-cr upper limit of 0.04 rayleighs. Similar limits apply to the Si III multiplet, the strongest 
members of which would be at 1108.36, 1109.97, and 1113.23 A. Phosphorus, which has about the 
same relative solar abundance as chlorine (Grevesse & Sauval 1998), also has strong transitions of 
its low ionization states, P II at 1154 A, P III at 918 and 1002 A, and P IV at 951 A. However, 
none of these multiplets are detected. The possible identification of an N III feature was mentioned 
above, but also likely to be present is N II, whose resonance multiplet is at 1085 A, but which is 
also not detected. 

The resonance transition of S V at 786 A has recently been reported by Steffi et al. (2003). 
There are only high-lying transitions of S V in the FUSE spectral range. An intercombination 
transition that connects to the ground state lies just beyond the long wavelength limit of the LiFlb 
detector at 1188 A. While there are also no transitions of either O II or O III connecting to their 
ground states there are several transitions between excited states of O II that could be observed. 
One such triplet, with its strongest component at 1083.13 A, has been identified in FUSE spectra 
of terrestrial airglow where its source is likely photoionization of neutral atomic oxygen (Feldman 
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et al. 2001b). The upper state of this transition lies 26.3 eV above the ground state and might be 
expected as we do observe S III A1126 and S IV A1099, which require 21.5 and 23.0 eV, respectively, 
for excitation (see Figure 2). The O II multiplet, however, like the N II resonance multiplet, falls 
in the gap between the two segments of each LiF detector and can be observed only with the SiC 
channels, whose effective areas at these wavelengths are a factor of five or more lower. There are 
also no collision strengths available for this O II transition. 

3.2. MDRS Spectra 

A portion of the higher resolution spectra corresponding to the top panel of Figure 1 is shown 
in Figure 5. The lower spectrum is from LiFla and contains some weak disk emission features in 
addition to the torus lines. The upper spectrum is from LiF2b and shows only torus lines but is 
considerably noisier due to the factor of two lower effective area of this channel relative to LiFla. 
In order to obtain the best S/N possible, all of the data from both dates of observation have been 
co-added for each of the channels. While the LiF2b spectra are useful in identifying torus features, 
only the LiFla spectra are used for fitting the line profiles. 

Of the six lines in the S III multiplet only the lines at 1015.496 and 1015.561 A are not 
resolved. The widths of the lines for this multiplet, the S IV multiplet at 1062-1073 A and the S II 
line at 1045.76 A are fit to gaussians and give linewidths of 0.14-0.15 A. The intrinsic linewidth 
for emission uniformly filling the MDRS aperture is evaluated from the terrestrial airglow triplet of 
O I at 1040 A, and found to be 0.085 A. The observed widths then correspond to the convolution 
of a gaussian of 7 — 8 x 10 5 K with the instrumental width for all of the three sulfur ions. This 
translates to an ion temperature of 60-70 eV. 

4. DISCUSSION 

4.1. Comparison with January 2000 Data 

In Paper I the spectrum was fit to an isothermal multi-level emission model using version 3 
of the CHIANTI database and software (Dere et al. 1997, 2001) that was updated with recent 
calculations of sulfur ion collision strengths that are now included in version 4 (Young et al. 2003). 
Normalizing to S III, the relative ion emissions were fit with an electron temperature, T e , of 80,000 K 
(6.9 eV), and relative abundances of S+ : S 2+ : S 3+ = 0.30 : 1.00 : 0.23. The present spectrum, in- 
cluding the spectral region below 995 A not observed in the January 2000 data, is again well fit with 
T e = 80,000 K, but with different relative abundances, 0.41:1.00:0.17. Comparing the ratios of 
the strongest lines in each observation, we find that relative to S 2+ , S + has increased by 37% while 
S 3+ has decreased by 26%. Further, comparing the relative ion intensities in the MDRS spectra 
from October 2001, we find that the relative abundances are about the same as in January 2001, 
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though this does not imply that they were stationary over the ten month interval separating the 
observations. 

Our results may be compared directly with those derived from a series of observations of the 
torus made by the Cassini ultraviolet spectrometer on 2001 January 6, 12, and 14, recently reported 
by Stem et al. (2003). From a composite spectrum from all three days spanning the wavelength 
interval from 650 to 1750 A, and using the current version of the CHIANTI database, they find 
an electron temperature of 4.6 eV and relative abundances of S + : S 2+ : S 3+ = 0.30 : 1.00 : 0.15 at 
6.3 They also simulate the suprathermal electron tail (^-distribution) by a superposition of 
three Maxwellians but find little difference in the derived abundances. As noted below in Section 4.2, 
the lower value of T e found by Steffi et al. is not surprising and is in fact in good agreement with 
the result from the reanalysis of the HUT data, also spanning a wider wavelength range than the 
FUSE spectrum, given in Paper I. The difference in derived relative S + abundance is likely due 
to the larger field-of-view of the FUSE aperture as Steffi et al. show its abundance to be steadily 
increasing inwards of 6.3 Rj. 



4.2. Electron Density and Temperature Diagnostics 

The far ultraviolet emissions observed by FUSE originate in the "warm torus", the region 
outward from ~5.8 Rj, where the plasma is characterized by an electron density, n e , ~2000 cm~ 3 , 
an electron temperature, T e , ~5 eV (1 eV = 1.16x 10 4 K), and an ion temperature, Tj, in the range of 
50-100 eV (c.f. Bagenal 1994). These parameters vary significantly with radial distance from Jupiter 
and also with height above the centrifugal plane as determined from both in situ measurements 
and remote long-slit spectroscopy, the latter primarily from ground-based observations of S II and 
S III (e.g., Thomas et al. 2001), which sample all of the torus. The remote observations exploit 
plasma diagnostics in the form of temperature or electron density sensitive spectral line ratios (Hall 
et al. 1994) and are independent of inferences based on the total flux in a given emission feature 
(Taylor et al. 1995). As noted above in Section 4.1, we use the CHIANTI database for the atomic 
data, primarily transition probabilities and collision strengths, for the evaluation of the requisite 
line ratios. 

There are, however, several caveats in using these diagnostics with the FUSE spectra. The 
relatively large field-of-view, 30" x 30", translates into a range of radial distance and height above 
the centrifugal plane of ±0.67 Rj over which these parameters are averaged as there is no spatial 
resolution within the aperture. This is in addition to the line-of-sight integration that is common 
to all remote observations. On the other hand, the derived parameters well characterize the regions 
where the emission is being produced. 

Another problem, which has been pointed out before, is that the use of theoretical, thermally 
averaged effective collision strengths assumes a Maxwellian electron velocity distribution whereas 
the torus plasma is known to have a suprathermal tail that could contribute significantly to the 
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excitation of the ultraviolet emissions. Sittler & Strobel (1987) found small departures from a 
Maxwellian distribution at Io's orbital distance with a suprathermal (T e ~1 keV) to thermal (T e 
~5 eV) electron density ratio of ~0.001. Further support for non-thermal electrons comes from the 
Galileo plasma science instrument (Frank & Paterson 2000). They found thermal core temperatures 
varying from 3 eV at ~6 i?j to 8 eV at 7.8 with departures from the core Maxwellian distribution 
for energies in excess of 30 eV. In an analysis of Ulysses spacecraft data for the hot torus at 7-9 
Meyer- Vernet et al. (1995) applied the theory of velocity filtration for non-Maxwellian plasmas 
given by Scudder (1992), to the centrifugal force potential constraining torus plasma. They used 
a generalized Lorentzian function, known as the "kappa" distribution, to represent a non-thermal 
distribution function that is quasi-Maxwellian at low energies but has a power-law decrease (with 
the exponent a function of re) at high energies. This distribution has been used to explain various 
features of both Voyager data (Moncuquet et al. 2002) and HST observations of Io's ultraviolet 
limb glow (Retherford et al. 2003). 

Thus, there is significant evidence in both in situ and remote sensing observations for non- 
Maxwellian plasma distributions in the Io torus in addition to the fact that Tj 3> T e . With the small 
re indices ~2-3.5 inferred by Meyer- Vernet et al. (1995), Moncuquet et al. (2002), and Retherford 
et al. (2003), the electron temperature increases from its value at the centrifugal equator, for 
example, 5 eV to ^7.5 and 15 eV at 1 plasma scale height above and below the equator for re = 3.5 
and 2, respectively. With the large aperture of FUSE, the derived effective electron temperature 
will depart significantly from its equatorial value. 

This effect is less important for the longer wavelength ultraviolet emissions, and evidence for 
this effect may be found in Paper I in the reanalysis of the torus spectra taken by the Hopkins 
Ultraviolet Telescope (HUT) during the Astro- 1 and -Emissions (Moos et al. 1991) using the current 
set of atomic data. For both of these data sets, the best fit to the emissions in the FUSE spectral 
range together with S II A1256, S III] A1729, and S IV A1406, gave a value of T e of ~50,000 K, similar 
to what was found earlier by Hall et al. (1994). Nevertheless, the derivation of an "equivalent" 
electron temperature is extremely useful for comparison with theoretical models of the electron 
energy distribution and with in situ measurements from Galileo. 

A final caveat rests with the atomic data, the uncertainties not being immediately quantifiable. 
There are a number of multiplets where the relative intensities within the multiplet as measured 
differ from the model synthetic spectrum and these lines are neither n e or T e dependent. There 
are also deviations from the predicted multiplet relative intensities, particularly for S II, but also 
for S III A1077.2 relative to the resonance multiplet. These deviations lie mainly in the details and 
may serve the purpose of being an "astrophysical laboratory" for fine-tuning the calculations. 

The strongest electron temperature dependence occurs for the multiplets of S III and S IV 
shown in Figure 2. As noted in Section 3.1.2, these multiplets require 21.5 and 23.0 eV, respectively, 
for excitation, while the resonance multiplets require only ~11 eV. Figure 6 shows the calculated 
T e dependence of three lines of S III and two lines of S IV, together with the observed ratios from 
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the LiFlb (crosses) and LiF2a (diamonds) channels. For all but one of these lines the two channels 
give nearly identical results. The discrepancy in the S III A1126.8 blend may be the result of a 
known calibration problem with the LiFlb channel (Sahnow et al. 2000). The data are fit by a T e 
ranging from 65,000 to 90,000 K. However, note that lines originating in the same multiplet give 
different values of T e , the difference well exceeding the statistical uncertainties in the measured 
line ratios. This limited set of data also does not provide any evidence that the S IV emissions are 
produced in a region of higher T e than those of S III, as Hall et al. (1994) had concluded from the 
HUT data. A complete reanalysis of the HUT spectra will be presented elsewhere. 

There are a limited number of electron density diagnostic intensity ratios available that arc 
relatively insensitive to T e . Figure 7 shows two for S II and S III, calculated with version 4 of 
CHIANTI for T e = 50,000 and 80,000 K, together with the observed ratios. All but one are 
insensitive to T e . The two S III ratios give widely differing values for n e , and all give densities 
in the range of 3000 - 5000 cm -3 , higher than the Voyager-based model value of 2000 cm -3 , but 
consistent with in situ Galileo PLS measurements of of 3600 ±400 cm~ 3 (Frank et al. 1996; Bagenal 
et al. 1997). The small slopes of the ratio curves make these calculated ratios extremely sensitive 
to small uncertainties in the atomic parameters and in fact there are significant differences in some 
of these curves between versions 3 and 4 of CHIANTI. 



4.3. Ion Temperatures 

The spectra presented above in Section 3.2 are the first ultraviolet spectra of the Io plasma 
torus with sufficient spectral resolution to allow the determination of ion temperatures from the 
measured line shapes. Thomas et al. (2001) have used ground-based high dispersion observations 
of S II and S III in the red and near infrared to determine T« as a function of radial distance from 
Jupiter. Unlike those measurements, the spectra described in Section 3.2 were obtained with the 
FUSE long-slit at ~1 Rj, oriented roughly normal to the planet's equatorial plane. This results in 
a factor of two loss in flux over placing the slit at the ansa (Sandel & Broadfoot 1982; Steffi et al. 
2003) but effectively integrates the emission over all radial distances. The derived ion temperature 
thus refers to the ions in their region of maximum emission, roughly 6.0-6.5 Rj. 

We found that the linewidths, which are determined mostly by the distribution of gyroveloci- 
ties, corresponded to perpendicular ion temperatures of 60-70 eV for all three sulfur ions. This is 
in excellent agreement with the values of Tj for S II and S III in the region from 6.0 to 6.5 i?j found 
by Thomas et al. (2001). However, it is in strong disagreement with the results from the Galileo 
plasma spectrometer which Crary et al. (1998) interpreted to give Tj = 10-24 eV over this range 
of radial distance. 

In addition to effects of non-Maxwellian plasma distributions along the magnetic field, the 
processes that add ion mass and energy to the torus plasma, ionization of neutrals that escape Io by 
electron impact and charge exchange, create a highly non-Maxwellian distribution of gyro- velocity. 



-10- 



Fresh ions are accelerated by Jupiter's co-rotational electric field from essentially Io's orbital velocity 
to co-rotation with Jupiter and acquire a gyro- velocity equal to the velocity difference. Because most 
ions are created at Io's orbital distance, the gyro-velocity distribution functions for newly created 
ions is approximately a delta function with only a narrow spread in velocity. Smith & Strobel 
(1985) demonstrated in a Fokker-Planck calculation of distribution functions for perpendicular (to 
the magnetic field) velocity that if the residence time of ions in the torus is short compared with 
the time to collisionally evolve into Maxwellian distributions, the ion distribution functions remain 
highly non-Maxwellian. For conditions at the time of Voyager- 1, they found that S II and O II 
departed significantly from Maxwellian distributions and that the average energies/temperatures 
were highest for O II and S II and coldest for O III and S IV with differences approaching a factor 
of 2. The latter ions had essentially Maxwellian distributions. Smith & Strobel fit Maxwellian 
distributions to the quasi-thermal cores of their calculated distribution functions and found common 
core temperatures for sulfur ions under Voyager 1 conditions of 60-74 eV. Only O II yielded a high 
core temperature of 105 eV. 

Our procedure of fitting gaussians to the observed line widths captures the quasi-thermal cores 
of the sulfur ion distribution functions. The data have insufficient S/N to permit detection of the 
non-thermal component, which in the Voyager-1 case, according to Smith & Strobel, 5% of S II 
were in a suprathermal tail with a distinct pickup ion signature. The ion temperatures of 60-70 eV 
inferred from the FUSE spectra are comparable to the Voyager-1 encounter ion temperatures 
derived by Bagenal (1994) from Voyager PLS data and with the theory of Smith & Strobel (1985). 



4.4. Relative Minor Ion Abundances 

As in Paper I, we can estimate the relative abundances of minor ions in the torus by comparing 
the excitation rates of these ion lines to that of the S III line at 1012.5 A. Repeating from that 
paper, the collisional excitation rate is given by: 



Ci = ^T(T e )e-*W 



T l/2 
J- e 

where T(T e ) is the electron temperature dependent effective collision strength, Ei is the energy 
of the upper state of the transition, and k is Boltzmann's constant. Values of T(T e ) are taken 
largely from the CHIANTI 4.0 database (Young et al. 2003) except for the recent calculations for 
CI III given by Ramsbottom et al. (2001) which is not in the current version of the database. 
These include the calculations of Tayal & Gupta (1999) for S III, Blum & Pradhan (1992) for C II, 
Berrington et al. (1989) for C III and Dufton & Kingston (1989) for Si III. We find that [Cl 2+ ]/[S 2+ ] 
= 0.021 ± 0.0035, or [Cl 2+ ]/(all sulfur ions) = 0.013 ± 0.0022. These are somewhat lower than the 
values of 0.031 and 0.020, respectively, given in Paper I. We still lack the atomic data for CI II 
to definitively show that Cl 2+ is the principal ionization stage of chlorine in the torus. Note that 
our derived abundance is comparable to the cosmic abundance ratio of [C1]/[S] = 0.015 given by 
Grevesse & Sauval (1998). As before, since sulfur can account for no more than half of the total 
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ion density in the torus this leads to the conclusion that chlorine ions are present at the level of 
1% or less of the total ion population. There appears to be a statistically significant decrease of 
~30% in the Cl 2+ abundance relative to that of the sulfur ions, between January 2000 and January 
2001. It also appears that the abundance of Cl + relative to Cl 2+ increased over this period. This 
variability may indicate a primarily volcanic source of CI in Io's atmosphere (Lellouch et al. 2003; 
Moses et al. 2002). 

From the detection of C III A977 discussed in Section 3.1.2, we find that [C 2+ ]/[S 2+ ] = 3.7 x 
10~ 4 . However, this should be regarded as an upper limit to the abundance of Iogenic carbon in 
that these ions could have a solar wind origin (Krimigis & Roelof 1983). From the C II and Si III 
upper limits, 3a relative column density upper limits are [C + ]/[S 2+ ] < 2.5 x 10~ 3 and [Si 2+ ]/[S 2+ ] 
< 4.1 x 10~ 3 . From Grevesse Sz Sauval the cosmic abundance ratios relative to S are [C]/[S] = 
15.5 and [Si]/[S] = 1.66, so assuming that C 2+ and Si 2+ are also the dominant ion stages, both 
carbon and silicon are severely depleted in the torus relative to solar abundances. For carbon the 
depletion is greater than four orders of magnitude and is likely due to the outgassing of volatiles 
from Io from continuous volcanic activity. For phosphorus an upper limit abundance estimate is 
not possible due to the lack of published collision strengths. 

5. CONCLUSION 

The Far Ultraviolet Spectroscopic Explorer has recorded the spectrum of the Io plasma torus in 
the range 905-1187 A at spectral resolutions of 0.26 and 0.085 A, the highest to date. Almost all of 
the emissions observed arise from ions of sulfur and the details of the spectra can be matched well 
by isothermal models using the CHIANTI database of atomic parameters and multi-level excitation 
models. These models also permit constraints to be set on the electron density and temperature 
in the region of the torus where the ultraviolet emissions are excited. In addition to sulfur, the 
resonance multiplets of CI II and CI III are clearly detected and we derive an overall chlorine ion 
abundance in the torus of slightly less than 1%, a decrease relative to that derived from observations 
made a year earlier (Feldman et al. 2001a). C III A977 is detected and a number of features near 
990 A remain unidentified. The source of the carbon ions is uncertain and may be the result of 
charge transfer of solar wind ions in the Jovian magnetosphere. There is no evidence for ions of 
nitrogen, silicon, or phosphorus, all which which have strong lines in the FUSE spectral range. 
From spectra obtained at 0.085 A resolution we are able to measure the linewidths of the strongest 
far ultraviolet emissions. These are found to be consistent with an ion temperature of 60-70 eV for 
each of the observed sulfur ions. 

We thank Peter Young for valuable discussions and insights into CHIANTI and Andrew Stem 
for making available preprints of the Cassini results. We thank the FUSE ground system personnel, 
particularly B. Roberts, T. Ake, A. Berman, B. Gawne, and J. Andersen, for their efforts in planning 
and executing these difficult moving target observations. This work is based on data obtained 
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Table 1: Ion Lines Identified in the FUSE Spectrum of the West Ansa, 14 January 2001 
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"Laboratory wavelength. 

6 Average brightness in 30" x 30" aperture. 

c Standard deviation in observed counts. 
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Table 2: Ion Lines Identified in the FUSE Spectrum of the West Ansa - continued 
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"Laboratory wavelength. 

6 Average brightness in 30" x 30" aperture. 

c Standard deviation in observed counts. 
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FIGURE CAPTIONS 

Fig. 1. — FUSE spectrum of the west ansa of the Io torus obtained on 2001 January 14 showing 
the regions containing the CI III (top) and CI II (bottom) multiplets. The total exposure time was 
9740 s. The gray curve is the data scaled by a factor of 10. 

Fig. 2. — Same as Figure 1 for the spectral region from 905 to 940 A showing the S II second 
resonance triplet. 

Fig. 3. — Same as Figure 1 for the spectral region from 1085 to 1135 A showing the highly excited 
S III and S IV multiplets. 

Fig. 4. — Same as Figure 1 for the spectral region from 970 to 1000 A showing spectra from three 
independent channels. Unidentified features that appear in at least two channels are denoted by 

Fig. 5.— FUSE spectrum from the MDRS aperture obtained on 2001 October 19 and 21. The 
same spectral region as the top panel of Figure 1 is shown from two separate channels. The LiF2b 
channel has about one-half the effective area as the LiFla channel. The total integration time was 
31,333 s. The LiFla spectrum also includes several H2 emissions from the Jovian disk. 

Fig. 6.— Calculated line ratios as a function of T e for S III AA1128.5, 1126.8, and 1122.4 (a, b, and 
c, respectively) relative to A1012.5 and S IV AA1099.5 and 1098.9 (d and e, respectively), relative to 
A1062.7. A plasma with spatially homogeneous, constant electron temperature defined by a single 
Maxwellian distribution and n e = 2000 cm~ 3 was assumed. The CHIANTI v. 4 software was used. 
The observed ratios are indicated by crosses (LiFlb) and diamonds (LiF2a). The error bars are 
1-cr in the observed counts. 

Fig. 7. — Calculated line ratios as a function of n e for S II and S III lines indicated in the figure 
for two values of T e , 50,000 K (solid line) and 80,000 K (dashed line). The CHIANTI v.4 software 
was used. The observed ratios are from the LiFla channel. The error bars are 1-cr in the observed 
counts. 
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Fig. 1. — FUSE spectrum of the west ansa of the Io torus obtained on 2001 January 14 showing 
the regions containing the CI III (top) and CI II (bottom) multiplets. The total exposure time was 
9740 s. The gray curve is the data scaled by a factor of 10. 
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Fig. 2. — Same as Figure 1 for the spectral region from 905 to 940 A showing the S II second 
resonance triplet. 
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Fig. 3. — Same as Figure 1 for the spectral region from 1085 to 1135 A showing the highly excited 
S III and S IV multiplets. 
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Fig. 4. — Same as Figure 1 for the spectral region from 970 to 1000 A showing spectra from three 
independent channels. Unidentified features that appear in at least two channels are denoted ?. 
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Fig. 5 — FUSE spectrum from the MDRS aperture obtained on 2001 October 19 and 21. The 
same spectral region as the top panel of Figure 1 is shown from two separate channels. The LiF2b 
channel has about one-half the effective area as the LiFla channel. The total integration time was 
31,333 s. The LiFla spectrum also includes several H2 emissions from the Jovian disk. 
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Fig. 6.— Calculated line ratios as a function of T e for S III AA1128.5, 1126.8, and 1122.4 (a, b, and 
c, respectively) relative to A1012.5 and S IV AA1099.5 and 1098.9 (d and e, respectively), relative to 
A1062.7. A plasma with spatially homogeneous, constant electron temperature defined by a single 
Maxwellian distribution and n e = 2000 cm" 3 was assumed. The CHIANTI v. 4 software was used. 
The observed ratios are indicated by crosses (LiFlb) and diamonds (LiF2a). The error bars are 
1-er in the observed counts. 
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Fig. 7. — Calculated line ratios as a function of n e for S II and S III lines indicated in the figure 
for two values of T e , 50,000 K (solid line) and 80,000 K (dashed line). The CHIANTI v.4 software 
was used. The observed ratios are from the LiFla channel. The error bars are l-cr in the observed 
counts. 



